Coulombic effects on magnetoconductivity oscillations induced by microwave 
excitation in multisubband two-dimensional electron systems 
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We develop a theory of magneto-oscillations in photoconductivity of multisubband two- 
dimensional electron systems which takes into account strong Coulomb interaction between elec- 
trons. In the presence of a magnetic field oriented perpendicular, internal electric fields of fluctu- 
ational origin cause fast drift velocities of electron orbit centers which affect probabilities of inter- 
subband scattering and the photoconductivity. For the electron system formed on the liquid helium 
surface, internal forces are shown to suppress the amplitude of magneto-oscillations, and change 
positions of magnetoconductivity minima which evolve in zero-resistance states for high radiation 
power. 

PACS numbers: 73.40.-c,73.20.-r,73.25.+i, 78.70.Gq 



I. INTRODUCTION 

Observation of microwave-resonance-induced magne- 
toconductivity oscillations in the multisubband two- 
dimensional (2D) electron system formed on the free 
surface of liquid heliuroii^ have attracted much inter- 
est. These oscillations appear to be very similar to 
l/S-periodic resistivity (Rxx) oscillations observed in 
ultrahigh-mobility GaAs/AlGaAs heterostructures sub- 
jected to a dc magnetic field and to strong microwave 
(MW) radiation of a quite arbitrary frequency uj > uj^"— 
(here Wc is the cyclotron frequency). At low tempera- 
tures T ~ 0.2 K, oscillations reported for surface elec- 
trons (SEs) on liquid helium have a peculiar shape 
which is reminiscent of a derivative of a peaky func- 
tion rather than a sum of simple maxima expected for 
usual scattering. With an increase in radiation power, 
the minima of the magnetooscillations evolve in zero- 
magnetoconductivity axx states^, which are very simi- 
lar to zero-resistance states (ZRS) observed in semicon- 
ductor 2D electron systems^i^, because, under a strong 
magnetic field, both axx and Rxx are proportional to the 
effective collision frequency ^{B). The only important 
difference of these two similar phenomena is that for SEs 
on liquid helium magneto-oscillations are observed only 
in the vicinity of a specific MW frequency u = uj2.i, rep- 
resenting the resonance frequency for excitation of the 
second surface subband. 

In semiconductor systems, ZRS are explained® as a 
consequence of the negative linear conductivity condi- 
tion Uxx < which appears for high radiation power. 
Due to instability of the system under this condition, it 
enters a nonlinear regime and develops a steady current 
state with (Jxx{io) — 0. Mechanisms of the negative lin- 
ear response conductivity of a 2D electron gas formed in 
semiconductor structures are based on photon-induced 
impurity scattering within the ground subband^*^. They 
were quite successful in explaining magneto-oscillations 
and ZRS induced by MW radiation in semiconductor sys- 
tems. Unfortunately, they cannot be applied for expla- 
nation of negative conductivity effects in the 2D electron 



system formed on the liquid helium surface, because the 
MW frequency w considered in these theories is not re- 
stricted by the condition u — lo2,i- 

The new mechanism of the negative linear response 
conductivity reported recently^i^^ is based on nonequi- 
librium filling of the second surface subband. 



N2 > Ni exp {-hu2,i/Te) 



(1) 



induced by MW excitation (here Ni is the number of 
electrons at the corresponding surface subband and Tg 
is the electron temperature). This condition provides a 
new channel for negative momentum relaxation due to 
usual quasi-elastic inter-subband scattering. This possi- 
bility can be seen already from the energy conservation 
for electron scattering from the excited subband (l — 2) 
to the ground subband {I = 1). 

Consider the electron energy spectrum in a magnetic 
filed applied perpendicular to the surface 



sin,x = A; + huc [n + 1/2) - eEi\X 



(2) 



(here A/ is the spectrum of SE states, n — 0,1,2... , X 
is the center coordinate of the cyclotron motion, and iSy 
is the dc electric field directed antiparallel to the a;-axis). 
Then, the energy conservation yields 



hojc {n - n') + nw2,i - eSii {X - X') = 0. 



(3) 



In this equation, hL02,i is not a photon quantum. It is the 
energy difference for electron excitation in the z-direction 
(huJi^i' = A/ — A/'). Since a;2.i is substantially higher than 
u!c, and ei^yLs ^ huic (here = hc/eB), scattering 
down the surface levels means scattering up the Landau 
levels n' — n = m* > and we have 



eE\\ {X - X') = hwc [ — - m 



(4) 



Thus, quasi-elastic inter-subband scattering from I — 2 
to ^ = 1 will be the scattering against the driving force 
{X' < X), if uj2.i/ujc — m* > 0, or when B is a bit 
lower than the level matching point uj2.i/^c — "m* ■ At 



2 



the same conditions, electron scattering up the surface 
levels is obviously the scattering along the driving force 
{X' > X). Therefore, for the appearance of the nega- 
tive conductivity correction, the additional condition of 
Eq. ([U is necessary. 

In semiconductor 2D electron systems^, resistance 
minima occur at uj/uJc = fn* + 1/4. For surface electrons 
on liquid helium of the areal density 12^ ~ 10^ cm~^, mag- 
netoconductivity minima reported in Ref. [3 are also lo- 
cated approximately near these "magic" numbers m* + 
1/4. Therefore, it is interesting to investigate if there 
is really such strong correlation in positions of Rxx and 
axx minima for these two different phenomena induced 
by MW excitation. In the single-electron theory of MW- 
induced magnetooscillationsi^, the conductivity minima 
are placed substantially closer to the level matching num- 
bers m*, because the Landau level broadening is very 
small for SEs on liquid helium. It is expected that 
larger distances of axx minima from the level matching 
points observed in the experiment will be explained by 
the many-electron effect. An additional interest in study- 
ing this problem is inspired by the observation of the res- 
onant photovoltaic effect which emerges at the minima 
of conductivity oscillationsii. 



The important conclusion of the single-electron treat- 
ment of MW-induced magneto-oscillations is that the 
negative conductivity terms are large for an electron sys- 
tem with extremely narrow Landau levels, when the col- 
lision broadening of Landau levels r„ is much smaller 
than temperature T. This condition is well realized for 
surface electrons on liquid helium. Still, under usual 
experimental conditions, this electron system is in the 
strong Coulomb coupling regime, which means that the 
average Coulomb interaction energy of an electron Uc is 
much larger than T and hujc- Therefore, a thorough in- 
vestigation of Coulombic effects on magnetoconductivity 
oscillations is required. 



In this work, we report the many-electron theory of 
magnetoconductivity oscillations induced by resonance 
MW excitation which takes into account strong Coulomb 
interaction between electrons. We use the fluctuational 
electric field conceplii^ which was generalizedi^ii^ to in- 
corporate the collision broadening of Landau levels. In 
the new theory, the Coulomb interaction between elec- 
trons affects strongly the line-shape and broadening of 
magnetoconductivity oscillations. As a result, positions 
of axx minima become dependent strongly on electron 
density rig and on the level matching number m*. This 
Coulombic effect appears to be very sensitive to electron 
temperature, therefore, it could be used for estimation of 
electron heating in experiments on MW-induced magne- 
tooscillations. 



II. DECAY RATE AND OCCUPANCY OF 
EXCITED SUBBANDS 

At low enough temperatures, depending on electron 
density, the 2D electron system formed on the free sur- 
face of liquid helium undergoes the Wigner solid transi- 
tion. This transition occurs when Uc /T ~ 137—. There- 
fore, there is a broad range 1 <C Uc/T < 137, where the 
electron system represents a highly correlated Coulomb 
liquid. These conditions arc quite usual for experiments 
on SEs in liquid helium. It is well established that al- 
ready at Uc/T > 9 the electron velocity autocorrelation 
time Tg-e is very high and close to the reciprocal of the 
harmonic oscillator frequency in a 2D triangular electron 
latticei^. This allows to simplify the transport theory, 
using the assumption that SEs are at equilibrium in the 
moving center-of-mass reference frame. In the semiclas- 
sical transport theory, this simplest many-electron effect 
increases the effective collision frequency of SEs by a nu- 
merical factor of the order of 2, which depends on the 
perpendicular electric field E±. In the presence of the 
quantizing magnetic field, another kind of many-electron 
effect becomes important which can substantially reduce 
the effective collision frequency as compared to the result 
of the single-electron treatment. 

The influence of Coulomb interaction on electron scat- 
tering under magnetic field can be well understood in 
terms of the quasi-uniform electric field of fluctuational 
origin, acting on each electron of a 2D electron liquid. 
It is assumed that, under the condition Uc 3> Te i an 
electron is displaced from an equilibrium position due to 
thermal fluctuations^^ similar to an electron in a triangu- 
lar Wigner lattice. When the magnetic length is much 
smaller than the electron spacing a, the restoring electric 

field of other electrons can be considered as a quasi- 
uniform field, causing a fast drift velocity of the elec- 
tron orbit center Uj' = c-B|''/i3 directed perpendicular 

to Ey. . The fast motion of cyclotron obit centers affects 
strongly the magnetoconductivity of electrons scattered 
by helium vapor atomsiS and capillary wave quanta (rip- 
plons)^^. Monte-Carlo simulations^® indicate that, in a 
wide range of the ratio Uc/T^ ^ 1, the distribution of the 
fluctuational electric fleld is close t o a Gaussian with the 

width parameter given by i^^f) = — '^\fT'en'J^ ■ 

It is remarkable that a strong quasi-uniform fluctua- 
tional field does not smear off Landau levels, because 
any uniform field can be eliminated by a proper choice of 
the reference frame. It is the frame where the cyclotron 
orbit center is at rest, and the electron spectrum coin- 
cides with the usual Landau spectrum. In the laboratory 
frame, the fluctuational fleld tilts the Landau levels simi- 
lar to Eq. We shall generalize the many-electron the- 
ory of quantum magnetotransport in 2D Coulomb liquids 
developed in Ref. [l3| and [3 by including in consideration 
inter-subband scattering. In particular, we shall consider 
the decay rate of the first excited subband, which deter- 
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mines subband occupancies under the MW resonance. 

The ZRS of SEs are observed at low temperatures T < 
0.2 K, where electrons are predominantly scattered by 
capillary wave quanta. The electron-ripplon interaction 
is usually described by the following Hamiltonian^"* 



Hint (R-e 



-T 

q 



Uq (z) = (qz) + eE^ - 



y y 



dVe 



(0) 



dZe 



(5) 



(6) 



where A is the surface area, Qq — ^Jhq/2pujq, u!q ~ 
y^a/pq'^^'^ is the ripplon spectrum, Rg = {ze,re}, hq 



q and 6q are the creation 



is the ripplon momentum, b 
and destruction operators, A = (e — 1) /4 (e + 1), e is 
the liquid helium dielectric constant, and Ki (x) is the 
modified Bessel function of the second kind. 

In the presence of the perpendicular electric field E± , 
the wave functions of the first two surface subbands can 
be approximated by 



ipi{z) — Aize 



ip2{z) A2Z 



1 



272,1 



(7) 

where Ai are normalization constants, parameters 7; are 
found by the variation, and 7;_;' = (7;+7;/)/2. In 
the magnetotransport theory, using the dimensionless 
parameter Xq = it is convenient to represent 

matrix elements of Eq. ^ as {Uq)^, 
where 



KVi^i. ixq)/Ll, 



where A 



i.i' 



[min(n,nO]! ^l„_„/|^_^ 
[max(n, n')]l 



V' 



I' I) 



(n.n' ) 



(10) 



mini n.n 



ix) 



L™(x) are the associated Laguerre polynomials, and Nq 
is the ripplon distribution function. 

Since Eq. (ITUl) is independent of X, one can find 
the average rate Di^i' using equilibrium electron dis- 
tribution over Landau levels: Z^^e~'^"^'^'' , where Z|| = 

^„ e~^"/"^=, and £„ is the Landau spectrum. For an en- 
semble of electrons, this procedure is equivalent to the as- 
sumption that electrons are in equilibrium in the moving 
center-of-mass reference frame, where E'| = 0. In order 
to include the collision broadening of Landau levels, we 
follow Ref. and represent the probability of scattering 
in terms of density-of-state functions of the initial and 
final states. Then, (5 (e — e„) of the density-of-state func- 
tions can be replaced by — ;^ImG;_„ (e). Here Gi^n (e) is 
the single-electron Green's function of the corresponding 
subband, whose imaginary part is broadened due to col- 
lisions with scatterers. This procedure is equivalent to 
the self-consistent Born approximation (SCBA)2£. 

Using the Gaussian shape of Landau level densities, 
inter-subband scattering rates can be represented as 



q 



xi.i' (q) sl°J iq,uji,u - q • Ud) 



(11) 



where TOo is the free electron mass. 



Xi,i' 



2ah?x, 



(12) 



Vi^i' (x) = xwi^i' 



l£±B_ 
A '■ 



(8) 



27ri/^fi ^ Jin 
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n.n' 



l,n:l' . 



F, 



Ll' 



Ll' 



1/2 



dv 

dz 



1/2 



(9) 



v(z) = —K/z -\- eE±z, and the functions wij> (y) repre- 
sent the matrix elements of 2W {qz). Useful expressions 
for the dimensionless functions wim (jj) are given in the 
Appendix. 

At first, consider the probability of inter-subband scat- 
tering of an electron moving in crossed magnetic B and 
electric Ey fields. The matrix elements of the interaction 
Hamiltonian yield X' — X ^ —qyL^, therefore, the term 
ei?|| {X — X') entering the energy conservation of Eq. ([3]) 
can be rewritten as ftq- u^, where is the drift velocity 
(ud = cE\\/B). Then, in the Born approximation, rates 
of inter-subband scattering can be found as 



yi,n^V,n' = ¥ E (2iV, + 1) {Uq) 



q 



X exp 



^ l.r 



r 



^ l.n 



(13) 



m — n' — n, Ti^n is the collision broadening of Landau 
levels of the corresponding subband, and 2r^,^.;, — 
^Ln + ^l.n'- In Eq. (HI]), we consider Nq ~ T/'hwq^ 1. 
For short-range scatterers, the collision broadening is in- 
dependent of n (r = h^j2uJcvlT^, here v is the collision 
frequency in the absence of the magnetic field). In the 
case of electron-ripplon interaction, the collision broad- 
ening of Landau levels depends on r^. At q • — > 0, 
Eq. (|lip coincides with the decay rate found previously 
in the absence of E||^*''-^^. 



For / 



V , the factor sf'i, {q,u}) represents the dynami- 



cal structure factor of a nondegenerate 2D system of non- 
interacting electrons. The definition of Eq. extends 
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the factor for I ^ V to describe inter-subband scatter- 
ing. Relationship (q, uj) = sl^J {q,uj — q - u^), which 
follows from Eq. pT|) . represents the Galilean invariance 
along the interface for the generalized factor, assuming 
electrons are in equilibrium in the center-of-mass refer- 
ence frame, where E'| = 0. 

The Eq. describes how a fast drift velocity affects 
the probability of inter-subband scattering. As noted 
above, the Coulomb liquid can be approximately con- 
sidered as an ensemble of independent electrons whose 
orbit centers move fast (u^; — >• u^'^) due to the fluc- 

(i) 

tuational electric field E} • In experiments on MW- 
induced magneto-oscillations of surface electrons in liquid 
helium, typical electron spacing a varies approximately 
from 1 • 10~^ cm to 3 • 10"** cm, which is much larger than 
the typical electron localization radius in the perpendic- 
ular direction ~ 10~^cm. Therefore, we can assume 
that the fluctuational electric field acting on electrons, 
occupying first few excited subbands, is the same as that 
acting on electrons of the ground subband. 

In the single-electron theory^", scattering probabilities 
and the momentum relaxation rate of a multisubband 2D 
electron system are found in terms of the generalized dy- 
namic factor S^^J (q, lu) using its equilibrium properties. 
Therefore, the main goal of the many-electron theory is 
to construct the Coulomb liquid version of this factor 
Slf,°^ iq,ui) applicable under conditions Uc/T ^ 1. Fol- 
lowing the ideas developed for intra-subband scattering, 
we expect that the many-electron dynamical factor of the 
multi-subband system of strongly interacting electrons, 
can be approximated as 



iq,u;)^(slf,iq,u;-q-Uf))^, (14) 



where {■■■) f denotes averaging over the fluctuational elec- 
tric field. For 1 = 1', this relationship was already 
proveni^ii^. For inter-subband scattering, the approxi- 
mation of Eq. (|14l) follows from the simple averaging of 
Eq. (|lip over the fluctuational electric field, where uj 

substitutes for u^. Since the maxima of S^^^, {q,i-o) as a 
function of frequency has a simple Gaussian shape, the 
model of Eq. provides larger broadening of these 



which agrees with the broadening of the dynamical struc- 
ture factor of the 2D Wigner solid under a strong mag- 
netic field--. 

When considering inter-subband scattering, we must 
separate scattering down and up the surface subbands. 
For electron scattering down [l > I'), at typical electron 
densities Ue ~ 10^ cm~^, the approximation of Eq. ([Ti)) is 
quite sufficient. For description of electron scattering up 
the surface subbands, an equilibrium form of S^'^'^^ {q, uj) 
should contain also small frequency shifts which provide 
the detailed balancing vi^i = Di^i' eyip {—Hujij> /Te). 

Correct frequency shifts of S^™,'^^ {q, lu) can be obtained 



from the detailed balancing, or using an accurate model 
based on properties of the 2D Wigner solid under a strong 
magnetic fieldMi^^. This approach yields 



Si I' (9. ^) = 



E 

n,n' 



i' + ^q^c 



(15) 



where 



llZ\n'i^)=^^P{-DLn-,l',n'{L0)}, (16) 



AT, 



AT, 



^ Ln:l' 



^ l,n 

8r2' 



and m = n' — n. Frequency shifts r^„/4Te and Xqr-^/4Te 
are small as compared to huc- It should be noted that the 
approximation of Eq. (ITi|) gives Eqs. (ITSI) and (IT51) with- 
out the small frequency shift XqT^/ATe- We restored it 
from the expression for the dynamical structure factor of 
the Wigner solid. This shift follows also from the detailed 
balancing, because it provides the important property of 
the equilibrium factor 



(17) 



which is equivalent to the detailed balancing. 

Thus, to describe inter-subband scattering rates vi^v 
in the multi-subband Coulomb liquid, one can use the 
general equation similar to Eq. (jlip . where the equilib- 
rium many-electron factor S^^f' {q,uji ii) of Eq. dTSl) sub- 
stitutes for sf'i, {q,uJi,i' — q • u^). The Coulomb broad- 
ening of the generalized factor ,Jx^Tc is not equivalent 
to the collision broadening of the SCBA theory, because 
it depends on the wave- vector argument of s['Z'^ {q,uj), 
which complicates evaluations. 

Consider the decay rate of the first excited surface sub- 
band V2^i- The usage of Eqs. ([T5|) and yields 



n 



V2- 



^ -D II n,m ^ ^ 



X J V{, (x) 




X exp 



ft^ (W2,l 



[L- {x)r , 

^2,n;l,n+m + C 



xTl/4hTe 



r2 

2,n;l,n+m 



-xTl 



dx, (18) 



where m is a positive integer, and the frequency shift 
r;^^/4re is disregarded as well as the small parame- 
ter r^^/8Tg. Here we still keep the Coulomb shift 
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FIG. 1. Decay rate of the first excited subband !^2->i vs the 
magnetic field B for T = 0.2 K and rie = lO'^cm"^: single- 
electron treatment (dotted); many-electron theory [Te = T 
(solid), Te = 0.5 K (dashed), and = IK (dash-dotted)]. 
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FIG. 2. The occupancy of the first excited subband vs the 
parameter a;2,i/i^c — m* for T = 0.2 K, = 2 ■ 10®cm~^: 
single-electron treatment (dotted), and many-electron theory 
[Q.R = 2.5 • 10^ s"^ (solid), and 0^ = 5 . 10^ s'^ (dashed)] . 



XqT^/AhTe because it can be important at high electron 
densities. 

For typical conditions of experimentsii^ performed on 
SEs over liquid '^He, the results of numerical evaluation 
of Eq. (IT5|) are shown in Fig. [1] The magnetic field 
range is chosen to be a vicinity of the level matching 
number m* = 4. One can see that already at a very 
low electron density = lO^cm"^, the many-electron 
curve calculated for Te = T = 0.2 K (solid line) sub- 
stantially broader than the result given by the single- 
electron theory (dotted line). Under these conditions 
Uc/T ~ 15. The broadening of the many-electron curve 
increases strongly with rig and the level matching num- 
ber TO* due to Fc c>c n'J^ j B^^''^ . Another important 
conclusion, which follows from this figure, is that the 
many-electron curve becomes even broader with heating 
of electrons because Fc oc s/Tl. It should be noted also 
that the single-electron decay rate has much weaker de- 
pendence on Tg which is not shown in Fig. [TJ 

The above noted Coulomb broadening of the decay rate 
i>2-n affects strongly subband occupancies n/ = Ni/N 
in the presence of MW radiation. For the two-subband 
model, which is valid for Tg ^ '2.K, we have 

n2 1 -|-exp(-fiw2,l/^e)i^2^.l/»'mw .,„^ 

— - 1 : U»J 

"1 1 -I- I/2^l/'rmw 

where r^w is the MW excitation rate. Under the reso- 
nance condition, rmw = f^|j/27mw, where 7inw is the half- 
width of the resonance, and is the Rabi frequency 
proportional to the amplitude of the MW field. Thus, 
the fiuctuational electric field leads to much broader 
variations of n2{B) near the level matching condition 
^2.1 l^c = Tn* than those given by the single-electron 



theory. 

For electron density rig = 2 . lO^cm^^ and Tg = 
T = 0.2 K, variations of the fractional occupancy n2 are 
shown in Fig. [2] When conducting numerical evalua- 
tions, we assume that the width of the MW resonance 
is due to inhomogeneous broadening 27mw = 0.3 GHa^, 
and the MW excitation rate rmw is independent of B. 
The many-electron lines calculated for to* = 3 and 7 at 
^Ir = 2.5 • 10^ s~^ (solid) are substantially broader than 
the corresponding single-electron lines (dotted). The in- 
crease of n2 with MW power is shown by many-electron 
lines (dashed) calculated at VIr — 5 • 10'' s""'^. Since 
the negative conductivity correction is proportional to 
n2—ni exp {—huj2,i/Te), it is clear that the many-electron 
effect should displace strongly conductivity minima. 



III. LINEAR MAGNETOCONDUCTIVITY 

According to the analysis given in the previous Section, 
probabilities of electron scattering in a multi-subband 
Coulomb liquid {Uc/T > 10) can be described by the 
many-electron dynamic factor s'f^^'^ {q, lo) of Eq. (jl5p in- 
stead of sf'ii {q, uj) of the single-electron theory. There- 
fore, using equilibrium properties of 5'}™°'' {q, oj), like that 
given in Eq. (jl7l) , and following the formal procedure de- 
scribed in details in Ref. 10 , one can obtain an equa- 
tion for the momentum relaxation rate similar to that 
found in the single-electron treatment. The effective col- 
lision frequency u is a sum of contributions from intra- 
subband (i^intra) and inter-subband (j/mter) scattering. 
The Vintci consists of normal terms (i'n), existing even at 
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ni = nil exp{—huJi^i'/Te), and anomalous terms i^a pro- 
portional to fii — riii exp {—fkoi ii/Tf.). Simple replacement 

s\';>{q,u:)^s\T^ (g,c.) yields 



^intra 



^E^' / ^,XiAQ)Slr^iq.O)dx,, (20) 

e , J 



UTS! 



l>l' 



POO 

/ XqXi.i' (q) Slf,''\q,UJlJ')dXq, (21) 

Jo 




bJ 



47r ^ 



X / a;5X/,i' (2;g)c?2;9, (22) 
Jo 



FIG. 3. The effective collision frequency due to inter-subband 
scattering I'lnter vs the parameter u)2,i/oJc. (sohd). Two con- 
tributions vn and va are shows by dotted and dashed lines 
respectively. 



were 

(23) 

and S'l^"^''^ {q,^) = ds'f^^^ (^i^) /dijj. The expression for 
i^A was represented in the form of Eq. (1221) . containing 
terms with / > V only, by interchanging the running in- 
dices. The similar procedure resulted in Eq. (PT|) employs 
the property of Eq. ([TT]). 
Using the relationship 



e 

(24) 

which follows from Eqs. (|15p - p7|) . the expression for 
(xq) can be rearranged to a more convenient form: 



-,/(mo) 



-,(mc) 



(x,) = E 



I Jji^n' i^'^q) I 



l.n\l' 



r{m.e) / \ 



;i(wi,r-m^,) + (r2„,-r2„)/8re 



^ l,n;l',n' ^ -^"J^ C 



(25) 

where m = n' — n. It is important that the Coulomb shift 
XgY^/Te does not enter the numerator of Eq. ([25|) . where 

(j'f' n' ~ -l^f n) /^^e '^^^ ^e disregarded for extremely nar- 
row Landau levels. Therefore, the vanishing point of each 
term entering va{B) is unaffected by the many-electron 
effect and practically coincides with the level matching 
condition W2,i = {n' — n) ujc- At the same time, posi- 
tions of minima and maxima of i^a (S) are very sensitive 



to electron density due to the broadening of the general- 
ized dynamic factor ^r^„.;' „/ + XqV^. Additionally, at 
high rie, they are slightly affected by the frequency shift 
XqTl/AT, entering I^'^''},.^. {loi,v). 

Typical behavior of Vintc-c and its components {va and 
j^n) as functions of the parameter uj2.i/^c is shown in 
Fig. 13] for sufficiently high electron density (ng = 1.5 • 

lO^cm"^, Vlfi = 5-10''s~^), when maxima of S'j™°'' ((7,w) 
begin to overlap due to the Coulomb broadening ,Jx^Tc- 
Oscillations of va{B) and i'j<i{B) are quite different. The 
normal contribution j^n is just a sum of simple peaks 
centered at the level matching conditions cj2,i — m*ujc. 
With an increase of the parameter lj2,i/'-^c, maxima over- 
lap strongly, and eventually the dependence j^n(-B) ac- 
quires a positive background. 

The sign-changing correction va originates from the 
derivative 5';^™°^ (9jw), and, in the vicinity of the level 
matching conditions, it is close to zero, if overlapping 
is small. At larger to*, the overlapping of maxima of 

ril'-n' i^i.i') results in a negative background of va{B), 
because the amplitude of va oscillations decreases with 
TO*, and, at a fixed _B, the negative contribution of 
smaller m* dominates. It is remarkable that the nega- 
tive background of va{B) is compensated by the positive 
background of vt^{B) and the entire contribution from 
inter-subband scattering fintcr — J^n + i^A (solid line) os- 
cillates nearly with the zero background. 

The conductivity equation, which follows from the 
momentum-balance equation methodic, has the usual 
Drude form, where, instead of the quasiclassical collision 
frequency vq independent of B, it is necessary to use the 
effective collision frequency v = i^intra + J^n + I'A- There- 
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(solid), and 0^ = 5 ■ lO'^s-^ (dashed)]. 



fore, under strong magnetic fields, Uxx is proportional to 
I'intra + t'N + ^^A- Without MW radiation, the influence 
of the strong Coulomb interaction on Gxx is illustrated 
in Fig. H for Tie = lO^cm-2 and T = 0.2K [Uc/T ~ 47; 
liquid '^He). Tlie result of the single-electron treatment 
based on the SCBA theory (dash-dotted line) was ob- 
tained from Eq. ([^0]) by setting ni — > 1 and Tc 0. 
When lowering the magnetic field, the many-electron line 
(dashed) moves from the SCBA line to the Drude line 
(i/ — >■ i/Q, dotted line). Transition to the Drude line oc- 
curs due to the Coulomb broadening of „/ (0) which 
triggers off scattering processes with |n' — n| > 0. The 
many-electron line of Fig. Stakes into account terms with 
\n' - n| < 9. 

In the presence of MW radiation, typical behavior of 
cTxx [B) is shown in Fig. U by the solid line calculated 
for = 3.5 • 10^ s~^, rie = 10^ cm^'^, and for a re- 
stricted field range covering the level- matching numbers 
3 < TO* < 9. This line illustrates evolution of oscil- 
lations and conductivity minima with lowering B. In 
contrast to the single-electron theory, where the ampli- 
tude of conductivity oscillations increases with to* (low- 
ering i?)2iii£, here the amplitude of oscillations decreases 
with TO* for sufficiently low magnetic fields, due to the 

Coulomb broadening of ^2'n°i-n-fm' In accordance 

with Fig. [2J conductivity minima of Fig. 3] are quite dis- 
tant from the level matching points. 

The dependence of positions of conductivity minima 
on the level matching number to* is illustrated in Fig. [S] 
for lower density = 10^ cm~^ and two MW powers. 
This figure indicates that the distance of conductivity ex- 



tremes of the many-electron lines (solid) from the point 
^2,1 l^c — Tn* increases fast with to*. For the comparison 
reason, the corresponding single-electron line (dotted) of 
TO* = 7 is also plotted. It is interesting to note that un- 
der the assumption — T , a,n increase in power (see 
dashed lines) reduces a little bit the distance of 
tremes from the level matching points. This behavior 
correlates with the behavior of corresponding n2 lines in 
Fig. [21 The electron heating induced by MW radiation 
acts in the opposite way, increasing the Coulomb broad- 
ening of //™p.„/ {^1,1') (the frequency shift XqV'^/ATe is 
independent of Tg) due to Tc oc -v/tT. Therefore, posi- 
tions and broadening of the conductivity minima in an 
experiment can be used for estimation of electron temper- 
ature realized in this system in the presence of resonant 
radiation. 

Positions of Uxx extremes depend also on rie, as il- 
lustrated in Fig. [5] representing results obtained for two 
electron densities: 5- 10^ cm~^, and 2- 10^ cm~^. One can 
see differences in variations of conductivity minima and 
maxima. Conductivity minima definitely become more 
distant from the level matching points with the increase 
in rig for all to* presented in the figure. Regarding con- 
ductivity maxima, with the increase in rig, they become 
more distant for three smallest level matching numbers 
TO* = 4, 5 and 6, while for larger numbers (to* = 8 and 9) 
they become even closer to the level matching points. For 
the transition number to* = 7, the position of the cor- 
responding maximum is practically unchanged. The to* 
becomes larger, if lower electron densities are compared. 
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IV. DISCUSSION AND CONCLUSIONS 

The expression for the generalized dynamic factor 
given in Eqs. ([T5|) and together with equations for 
the effective collision frequency [Eqs. (|20 |) - ([22l) ] allows to 
describe photoconductivity of SEs on liquid helium in a 
wide range of the parameter Uc/T up to the Wigner solid 
transition. The main conclusion of the many-electron 
theory reported in this work is that the unusual MW- 
induced magneto-oscillations of <7xx and ZRS can exist 
in multisubband 2D electron systems even under strong 
Coulomb coupling conditions {Uc ^ T) relevant to ex- 
periments^. 

Coulomb forces acting between electrons affect 
strongly the broadening of magneto-oscillations and po- 
sitions of conductivity extremes, which become depen- 
dent on electron density, the level matching number m* 
and the electron temperature. Therefore, for SEs on liq- 
uid helium, in the many-electron theory, there are no 
fixed "magic" numbers describing positions of conduc- 
tivity minima on the a;2,i/wc-axis. The proximity of ex- 
perimentally observed minima to uj/uJc ~ m* -I- 1/4, most 
likely, is accidental (valid only for a specific electron den- 
sity). A substantial dependence of the distance of axx 
minima from the point ui /ujc = m* on the level matching 
number m*, which can be seen in experimental data of 
Ref. [3, agrees qualitatively with this conclusion. 

Numerical calculations of SE magnetoconductivity 
shown in the previous Section were performed under the 
assumption that electrons are not substantially heated 
by the MW (Tg ~ T) . This assumption allows the many- 
electron theory to reproduce remarkable shapes of a^x 
oscillations observed for SEs on liquid helium. The inclu- 
sion of Coulomb interaction in the theory affects crucially 



the dependence of the amplitude of magnetooscillations 
on the magnetic field strength. In the many-electron 
theory, the amplitude of oscillations eventually begins 
decreasing with lowering B, which agrees with experi- 
mental observations. Even though the Coulombic effect 
suppresses amplitudes of magneto-oscillations, under ex- 
perimental conditions of Ref. d, the linear theory results 
in axx < at certain ranges of magnetic and MW fields. 

The comparison of theoretical and experimental 
o'xxiB) curves indicates that typically the real heating 
of SEs is not strong and most likely ~ T. Since 
i^intra ~ l/T'e (Ff „, + x^r^)!/^ a stroug iucreasc in 
electron temperature (T^ ^ T) would lead to deep drops 
of axx at the both sides of the level matching conditions, 
which contradicts to observations. In the single-electron 
theory, with an increase in electron temperature the sign- 
changing correction va becomes more efficient than i^intra 
and i/N, because r;_„ is independent of Tg. In the many- 
electron theory, the Coulomb correction to the broaden- 
ing of sl™,'^^ {q,u!ij') increases with electron temperature, 

y^x^Tc oc Te , and the relative increase in efficiency of 
z/A as compared to I'intra and z^n approaches a saturation, 
when y/x^Tc S> Ti^n-i',n'- The electron heating affects 
positions of axx extremes, which also can be used for es- 
timation of Tg. The many-electron theory of cold SEs re- 
sults in minima distances from the level matching points 
which are of the same order, as those observed experi- 
mentally. This also agrees with the estimation ^ T 
given above. 

Regarding negative conductivity values obtained in the 
linear theory, regions with axx < are surely unsta- 
ble, because any electron density fiuctuation produces 
an electric field which increases the fluctuation. For the 
Coulomb liquid formed on the liquid helium surface, even 
the fluctuational field E^*^ can be a cause for instability, 
because excited SEs will be scattered against the restor- 
ing force — eE^ . The edge of the 2D electron system on 
liquid helium is usually fixed by application of a dc con- 
fining field of Corbino and guard electrodes which also 
can be a cause for instability of the system at axx < 0. 

A stable state can be achieved under the nonlinear 
transport regime by forming a strong steady current jo, 
of which axx{jo) = 0. In the 2D Coulomb liquid on the 
surface of liquid helium, formation of current domains is 
unlikely because of strong electron correlations. The SEs 
usually have no source and drain electrodes. Therefore, 
a steady current can be formed only by electrons circling 
the center of the electron pool. Since it is impossible to 
create a strong current density in the center, electrons 
will move to the edges of the electron liquid depleting 
the center. As a result, a nonuniform electron density 
distribution along the surface will be formed to provide 
radial electric field and a circling current, strong enough 
to make axx = 0. At a fixed magnetic field, a change 
in electron density additionally helps the system to leave 
the unstable regime, due to the many-electron effect de- 
scribed here. 
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Thus, under the condition a^x < of the hnear trans- 
port regime, electron moving against the confining force 
of Corbino and guard electrodes is the most likely way 
to reach a stable state. This scenario agrees with re- 
cent observations of the resonant photovoltaic effect for 
SEs on liquid heliun*ii, where MW excitation causes 
a strong displacement of surface electrons towards the 
edges (against the confining force) , if the magnetic field 
is close to the minima of the conductivity oscillations. 
Therefore, experimental observation of the strong density 
dependence of positions of conductivity minima consis- 
tent with the results of the many-electron theory would 
be an additional evidence for the negative conductivity 
as the origin of the resonant photovoltaic effect. 



where Ai,i = 1, and 



1 72,1 
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A2,l (x) = ( 1 - -X ) , A2,2 (x) 



The integrals of the right part of Eq. (|26l) can be found 
analytically, and we have w\^\ (y) — wi{y), 



W2,i (y) 
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Here we use the following notations 
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According to definitions given in Eqs. ©-([S]), for the 
two-subband model [l — 1,2), the functions wi^i' (y), de- 
scribing electron-ripplon coupling, can be written as 



wiM (y) 



47?,' 



Xij' (x) dx, 
(26) 



W2iy) 



^-y ^ 

(l-y) 
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wsiy) = 7z r^^ 1 ^In 



a-yr (l-y) 
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Vy 



Vy 



The analytical presentation of wiji (y) given above is very 
useful for numerical evaluations of a^x ■ 
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